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This report presents the computational results of a comparison of maser materials for 
a 32-GHz maser amplifier. The search for a better maser material is prompted by the rela- 
tively large amount of pump power required to sustain a population inversion in ruby at 
frequencies on the order of 30 GHz and above. The general requirements of a maser 
material and the specific problems with ruby are outlined. The spin Hamiltonian is used 
to calculate energy levels and transition probabilities for ruby and twelve other materials . 
A table is compiled of several attractive operating points for each of the materials ana- 
lyzed. All the materials analyzed possess operating points that could be superior to ruby . 
To complete the evaluation of the materials , measurements of inversion ratio and pump 
power requirements must be made in the future. 


I. Introduction 

This report describes the results of a theoretical evaluation 
of several paramagnetic materials being considered for use in a 
32-GHz maser amplifier. Previously, ruby has been very suc- 
cessfully employed in 2.3- and 8.4-GHz masers [1] , and in an 
18- to 26-GHz tunable maser [2] . However, due to a mono- 
tonically decreasing inversion ratio above 12 GHz for ruby (for 
push-pull pumping) ruby becomes less favorable at higher fre- 
quencies. The inversion ratio is defined as the ratio of the 
inverted-spin population difference with the pump on to the 
thermal-equilibrium population difference with the pump off, 
and is determined experimentally by the ratio of gain (dB) 
with the pump on, to absorption (dB) with the pump off. 
Moore and Neff [3] and Shell (private communication) mea- 
sured an inversion ratio of 1.1 at 32 GHz. In a recently com- 
pleted reflected-wave maser [5] , the inversion ratio was esti- 
mated to be 0.7 to 0.8. This low inversion ratio was a contrib- 


uting factor to the reduced gain-bandwidth of the maser com- 
pared to a similar 22-GHz maser [31] in which the inversion 
ratio was at least 1.6. Theoretical calculations of spin-lattice 
relaxation rates [29] suggest that the low inversion ratio 
(/ « 1) is inherent in ruby at this operating point (‘operat- 
ing point’ refers to a given dc field strength, crystal orienta- 
tion, and pumping scheme). Hence, other paramagnetic 
materials, as well as other operating points of ruby, are inves- 
tigated as a first step toward finding the best maser material 
at 32 GHz. 

Before proceeding with the materials evaluation, two points 
should be made. First, the inversion ratio in the 32-GHz re- 
flected-wave maser (RWM) could probably be improved by 
either using more pump power or by using a maser structure 
with a higher 0 at the pump frequency. In the former case, 
heating of the maser due to microwave losses will degrade the 
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gain and noise temperature performance, but to what extent 
is unknown. In the latter case, the resonant structure would 
significantly reduce the tunability of the RWM, but would 
improve the pump power coupling over a still useful instanta- 
neous bandwidth. However, even with these engineering 
improvements, the inversion ratio would still only approach 
unity. 

The second point is that the inversion ratio, even though 
central to the evaluation of a maser material, cannot be 
accurately addressed by the methods presented here. It is 
believed that the low inversion ratio of ruby at the 32-GHz 
RWM operating point is a result of an unfavorable set of spin- 
lattice relaxation rates. These relaxation rates and the corre- 
sponding inversion ratios have been calculated, but, due to the 
complexity of the calculations, will be presented in a future 
report. 

In Section II the requirements of a maser material are dis- 
cussed, including the difficulties with ruby. A list is provided 
of the materials evaluated. In Section III the use of the spin 
Hamiltonian to calculate energy levels and transition proba- 
bilities is outlined. In addition, limits on the inversion ratio 
and a material figure-of-merit are discussed and the method of 
computation is reviewed. In Section IV a table of several 
promising operating points is presented for each of the mate- 
rials analyzed. The conclusions are presented in Section V. 
The Appendix contains a table of measured relaxation times 
obtained from the literature for several materials of interest. 


II. Maser Material Considerations 

The choice of a maser material is the single most important 
factor in maser design. The maser material consists of a non- 
magnetic crystalline lattice lightly doped (0.01-0.1 percent) 
with paramagnetic ions. Detailed discussions of suitable para- 
magnetic ions and host crystals are given elsewhere [6] , [7] . 
Here, only an outline of the necessary and desirable material 
properties is given. 

The most common paramagnetic ions are transition metals 
and rare earths, because of their unfilled 3 d and 4/ electron 
shells. To operate in CW mode, the ion should possess an orbi- 
tal ground state with three or more spin levels. This eliminates 
most of the rare earths. The additional requirement of a negli- 
gible nuclear magnetic moment (a source of inhomogenous 
broadening) reduces the possible ions to Cr 3+ , Fe 3+ , Ni 2+ , and 
Gd 3 *. (Actually, it is not clear that such broadening would 
adversely affect a maser with an inhomogenous applied field.) 
Only Cr and Fe are considered here, as they are by far the 
most common choices of active ions. The electronic configura- 
tions of Cr 3+ and Fe 3+ are 3d 3 and 3 d 5 , respectively. The cor- 


responding free-ion ground states are 4 F 3/2 and 6 S 3/2 . Hence, 
Cr 3+ has four spin levels and Fe 3+ has six. 

The host crystal must be non-magnetic, non-metallic, and 
available in large single crystals of a high degree of perfection. 
The material should have a sufficiently high thermal conduc- 
tivity and small loss tangent at liquid helium temperatures and 
microwave frequencies to minimize heating of the lattice. To 
facilitate the microwave engineering, the crystal should possess 
a relatively isotropic and temperature-independent dielectric 
constant. Finally, the material should be machinable and 
chemically stable, and be able to withstand thermal cycling 
between liquid helium and room temperatures. 

The active ion substitutes for one of the metal ions in the 
crystal lattice. The local crystal electric field seen by the 
ion splits the highly degenerate orbital ground state into de- 
generate pairs (assuming the number of spin levels is even). 
This splitting of spin levels due to the crystalline electric field 
is termed the zero-field splitting (ZFS). The ZFS must be large 
enough to permit pump-induced transitions between non- 
adjacent spin levels. As a rule-of-thumb, the ZFS should be of 
the same order of magnitude as the signal frequency. 

Another very important material parameter is the spin- 
lattice relaxation time, which describes how long spins remain 
in an excited state before returning to thermal equilibrium 
with the lattice. At microwave frequencies, spin-lattice inter- 
action is the dominant spin relaxation mechanism. For rela- 
tively low concentrations of paramagnetic ions (<0.05 percent) 
and at liquid helium temperatures, the most significant spin- 
lattice interaction is thought to be the Kronig-Van Vleck 
mechanism [8] , in which lattice vibrations induce transitions 
between spin states, and spin-spin interactions are neglected. 
The spin-lattice relaxation times of the various transitions 
must be long enough (> msec) to permit saturation of the 
pumped levels with a reasonable amount of pump power. 
Impurities, ion clustering, and dislocations (all of which are a 
function of the crystal growth procedure) can shorten relaxa- 
tion times, so pure, defect-free crystals are preferable. 

Another material parameter is the (unbroadened) linewidth 
Af L of the material. At liquid helium temperatures, Af L is 
determined primarily by spin-spin interactions and is inversely 
proportional to the spin-spin relaxation time, which is the 
average length of time between random dephasing “collisions” 
of neighboring spins [9] . The linewidth is usually within the 
range of 10-100 MHz for solid-state maser materials. For linear 
stagger-tuned masers of bandwidth » Af L , Af L can be shown 
to have no first-order influence on the gain-bandwidth proper- 
ties of the maser (see Section III). However, if the taper is 
along the length of the material, a material with a smaller Af L 
may exhibit a larger noise temperature at one end of its band- 
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pass [9] . At present, it is not understood how Af L impacts 
pump power requirements. 

Finally, many materials possess two or more magnetically 
inequivalent sites (i.e., sites having different spectra) for the 
active ions to occupy, thus decreasing the density of useful 
ions. (Gain in dB is proportional to spin density.) In most 
materials, certain orientations exist for which these sites be- 
come equivalent. For such materials, only these orientations 
will be analyzed. 

As mentioned in Section I, difficulties were experienced 
using ruby at 32 GHz. A relatively low inversion ratio was ob- 
tained for the orientation employed in the 32-GHz RWM. The 
pump transitions in ruby at this orientation are quite weak, 
making it difficult to saturate the levels. This is a result of the 
small zero-field splitting (ZFS) of ruby: ZFS = 1 1 .4 GHz, which 
is only about one-third of the signal frequency. To obtain 
sufficient separation between spin levels for amplification at 
32 GHz, a relatively large magnetic field (11.8 kG) must be 
applied. This field becomes the dominant influence on the 
spins, far exceeding the effects of the local crystal field. Under 
such conditions, the spins in the lattice assume nearly pure- 
spin characteristics, as if the Cr ions existed freely in the 
magnetic field. The selection rules of quantum mechanics 
allow transitions only between adjacent pure-spin states [10] , 
thus leading to a small stimulated transition probability for 
the pump transitions in ruby. A material with a larger ZFS will 
in general have stronger pump transitions and will therefore 
better absorb the pump power, all other factors remaining 
the same. 

The materials analyzed in this work are listed below. Details 
of the crystal structure, orientation of magnetic axes, spin 
Hamiltonian, and site equivalence are given in the references. 


Ruby (Al 2 0 3 :Cr) 

[11] 

Emerald (Be 3 Al 2 Si 6 0 18 :Cr) 

[12] 

Spinel:Cr (MgAl 2 0 4 :Cr) 

[13] 

YAG:Cr (Y 2 Al s 0 12 :Cr) 

[14] 

YGG:Cr (Y 2 Ga 5 0 12 :Cr) 

[14] 

RutilerCr (Ti0 2 :Cr) 

[15] 

Zinc Tungstate :Cr (ZnW0 4 :Cr) 

[16] 

Andalusite:Cr (Al 2 SiO s :Cr) 

[17] 

Yttrium Oxide :Cr (Y 2 0 3 :Cr) 

[18] 

Rutile.Fe (Ti0 2 :Fe) 

[19] 

Zinc Tungstate:Fe (ZnW0 4 :Fe) 

[20] 


Andalusite:Fe (Al 2 Si0 5 :Fe) [21] 

Sapphire:Fe (A1 2 0 3 :Fe) [22] 

III. Evaluation of the Spin Hamiltonian and 
Associated Parameters 

The spin Hamiltonian H s describes the interaction of the 
electron spin of the paramagnetic ion with the local crystal 
field and with the applied magnetic field. Evaluation of H s 
allows calculation of the energy levels and transition probabili- 
ties for the spin system of a singlet orbital ground state ion in 
a radiation field. Detailed discussions of the spin Hamiltonians 
may be found in the literature [8] , [23] . Berwin [24] gives a 
detailed derivation of the spin Hamiltonian based on the for- 
mulation of Bleaney and Stevens [23] for ruby. In [25] , 
Berwin discusses the spin Hamiltonians for several of the 
materials listed in Section II. 

For most materials of interest, the assumption of an “inter- 
mediate” crystal field is made, meaning that the interaction 
energy of the crystal field with the ion falls between the 
Coulomb and spin-orbit interaction terms. In deriving H s , the 
spin-orbit and Zeeman terms are treated together as a pertur- 
bation on the singlet orbital ground state. The absence of 
orbital degeneracy is sufficient for quenching of the orbital 
angular momentum L [8] . That is, to first-order,! is equal to 
zero, so the ion behaves in a pure-spin-like manner. Up to a 
second-order perturbation, the spin -orbit coupling admixes the 
singlet ground state with higher-lying orbitals, restoring some 
of the orbital angular momentum. This second-order effect is 
the source of the ZFS. 

For ruby, the spin Hamiltonian has the form [11] 

H s = gpB-S + DSl ( 1 ) 

where the spectroscopic splitting-factor is g ^ 2 (g = 2.0023 
for pure spin), 0 is the Bohr magneton, B is the applied dc 
magnetic field, and 

S = xS x +yS y + zS t (2) 

is a vector of spin operators. The spin operators describe the 
observable properties of the paramagnetic ion spin states and 
may be conveniently written in matrix form [10] . The Carte- 
sian directions are along the principal axes of the magnetic 
complexes. These axes, usually expressed in terms of the crys- 
tallographic axes, are used to describe the symmetry^ of the 
magnetic resonance spectrum. The orientation of B is ex- 
pressed by the usual azimuthal and polar angles, 0 and 6. The 
constant D determines the ZFS and reflects, in principle, the 
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extent of admixing with higher lying orbitals. This spin Hamil- 
tonian exhibits axial symmetry about the magnetic z-axis and 
has a ZFS = 2\D\. 

For magnetic complexes of lower symmetry, additional 
spin operator terms may be needed to accurately specify the 
resonance spectrum; e.g., - S 2 y ) and (S* + S'* + sf ), the 

so-called orthorhombic and cubic terms [26] . The form of the 
required spin operator terms can sometimes be determined 
from crystal field theory through the use of equivalent oper- 
ators, as discussed in [24] . However, the crystal field approach 
usually assumes ionic bonding and neglects covalency effects. 
The fact that the coefficients of H s , and quite often the form 
of H s , must be determined experimentally (by fitting to EPR 
data) is quite likely a result of this assumption. 

Given an H s , one can solve for the spin energy levels and 
eigenstates by solving: 


The vector components are defined by the relation: 


a.. = 4-(a. . x + B. .y + 7. . z) 
u 2 v */ *r } 


( 7 ) 


where a, /3, and 7 are in general complex numbers. Noting 
Eq. (4), a, 0, and 7 are easily evaluated for any number of spin 
levels; in [24] the results for a 4-level spin system are given. 


As a means of comparing the relative strength of transitions, 
Siegman [9] has defined the quantity 


H* • 00 * • 


a - 


i// 1 1 2 


( 8 ) 


The stimulated transition probability can then be written as: 

w /,- = b 2g(f) 2 (9) 


H s | * ( > = E. | V/,.> / = 1,2, ... , 25+ 1 

( 3 ) 

where | \jj.) is usually written as a linear combination of pure- 
spin states, 


ItfO = a.\S) + fc,|S-l>+ ... +r.|-£> (4) 

and where E t is the energy of the ith level, and S is the spin of 
the ion. The labeling of the pure-spin states is identical to the 
labeling of states for a large applied field. Equation (3) is most 
easily solved by expressing it in matrix form and using the 
usual matrix methods to solve for the eigenvalues and eigen- 
vectors of H s . Note that H s is Hermitian. 

Knowing the eigenstates, one may then calculate the rate of 
stimulated transitions due to an RF magnetic field H x . Apply- 
ing Fermi’s golden rule, the probability of a transition between 
states / and j is [9] : 

w , 7 =±y 2 g(J)\(t i \g&S-H 1 \t / >\ 2 ( 5 ) 


The absolute maximum value of a 2 is obtained by using RF 
fields that are polarized such that H 1 is paralled to o. This 
value of a 2 is given by the trace of a a* [9] , which has the 
value: 


(o 2 )m BX = j(I<*I 2 + I0I 2 + IyI 2 ) (10) 


This expression will be used to compare transition strengths. 
This formulation does not consider whether or not the pre- 
scribed polarization of H x is achievable in a given microwave 
circuit. 


The gain of the maser is obviously of central importance 
and involves several important material parameters. The gain 
in dB of an unbroadened RWM or TWM [9] is 


G 


dB 



( 11 ) 


where s is the slowing factor, 2 is the maser length in free-space 
wavelengths, and Q m is the magnetic Q of the maser material. 
The magnetic Q is defined as 


where 7 = gfaj ft, and g(f) is the line shape (as a function of 
frequency) for the transition. The term sandwiched in the 
matrix element is the magnetic dipole interaction energy. 
Since non-spin operators may be removed from the matrix 
element in Eq. (5), it is convenient to define the vector 
quantity 

Ojj = IS I ( 6 ) 


_ energy stored in material 

^ m energy emitted per cycle by material 

( 12 ) 

Assuming « 1 (|f * 1/3 for /= 32 GHz, T= 4.2 K), the 
reciprocal of Q m may be expressed as [9] 

1 7 2 h hf Io 2 r) N 

Qm Wo kT no - of l eve l s ’ 
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where / is the inversion ratio, r? is the filling factor, and N is 
the spin density. The filling factor accounts for the fraction of 
RF field in the material and the degree to which the field is 
optimally polarized; its value varies between 0 and 1 . The spin 
density is determined by the concentration of paramagnetic 
ions; for 0.05 percent Cr concentration in ruby , jV= 2.35 X 1 0 19 
spins/cm 3 . From Eqs. (11) and (13), we may express G dB in 
terms of material parameters as: 


G 


dB 


Io 2 N 

A f L • no. of levels 


The inversion ratio is defined as 


(14) 


I = 


An 




AN.. 

) i 


05 ) 


where A Nj / = Nj - N i is the thermal-equilibrium population 
density-difference, and An (J is the population-density differ- 
ence under pumped conditions. To determine / for the signal 
transition of a multi-level spin system, one must solve a set of 
rate equations that accounts for both stimulated transitions 
and spin relaxation [9] . These rate equations may be simpli- 
fied by assuming steady-state conditions and saturated pump 
levels, and by neglecting the influence of the incoming signal. 
Since the relaxation rates are not known (this issue will be 
addressed in the aforementioned future report), the rate 
equations cannot be solved for the actual inversion ratio. 
Instead, assumptions are made about the relaxation rates and 
the corresponding / is determined. In one case, all relaxation 
rates are assumed to be equal and the inversion ratio is defined 
as l equaV An upper limit can be put on I by assuming an 
optimum set of relaxation rates exist [9] . This is defined as 
I opt . Note that the actual inversion ratio may be less than 
both I op t and I equai , but for pumping schemes employing two 
pumps, the actual inversion ratio often lies between I opt and 
I equaV Expressions for I opt and I equal are easily derived from 
the rate equations but will not be given here. 

A maser material “figure-of-merit,” indicating the gain- 
bandwidth potential of the material, was proposed [27] to be 
AfJQm • From Eq. (13), and considering only material param- 
eters, we find that 


A 4 a io 2 N 
O no. of levels 

m 


( 16 ) 


The computer analysis was performed on the JPL UNIVAC 
(F-system) with a modified version of an existing Fortran code 
[25] originally written for the analysis of ruby. The JPL 
Fortran subprogram HERMQR 1 was used to compute the 
eigenvalues and eigenvectors of H s . The existing code was 
modified to include the spin Hamiltonians of the other mate- 
rials. Other small modifications were also made. 

The inputs to the program are as follows: 

(1) material 

(2) range of 6 and 0-increment 

(3) a single value of 0 

(4) range of B and 5-increment 

(5) signal frequency window f L 0 , f HJ 

(6) minimum allowed value for max (a, 0, 7) for signal 
transition, a min 

The code is run for a given material at a given 0. Typical values 
are 6 = 0 to 90 deg, Ad - 10 deg, 0 deg < 0 < 90 deg, and 
B = 0 to 15 kG, AB = 0.5 kG. For a 32-GHz signal frequency, 
the window was usually f LO = 31 GHz and f HI = 33 GHz. For 
transitions falling within this range, o min = 1.0 was chosen. 
The signal frequency window and o min are used to pre-select 
operating points, thus decreasing computer output. 

The program outputs are the following: 

(1) energy levels and eigenstates 

(2) transition probabilities (both a, 0, 7 and a 2 ) 

(3) local values for— ~ and 

V ' AS A f s 

(4) 1 op V 1 equal* and figure-of-merit 

Because of the large number of possible pumping schemes 
(especially for the 6-level systems), the following guidelines 
were employed in choosing schemes: 

(1) signal transition is between adjacent levels 

(2) use two pumps, when possible 

(3) pumps can skip one level at most 

These guidelines limit the pumping schemes to the usual ones 
employed. Many other schemes are possible [9] . 


Hence, the figure-of-merit optimizes the product G dB • A f L . 
For the case of a linear stagger-tuned maser with bandwidth 
Af A f L , one can show that the quantity on the right-hand 

side of Eq. (16) optimizes the product G dB • Af. 


l JPL Fortran V Subprogram Directory , Fifth Edition, JPL Publication 
D-829 (internal document), Jet Propulsion Laboratory, Pasadena, 
California, July 1982. 
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IV. Results and Discussion 

Table 1 compares several of the more promising operating 
points of each of the materials analyzed, starting with ruby. 
The materials are arranged in order of increasing ZFS; for 
materials with S = 5/2, the ZFS of the lower degenerate states 
is used. Across the top of the table is the material name and 
operating point number. The first row of the table shows the 
paramagnetic ion used. The second row gives the ZFS. Note 
that all the materials have a ZFS larger than ruby. The third 
row indicates the number of magnetically non-equivalent ionic 
sites in the lattice; a “1” means that all sites are equivalent. 
The fourth row gives the orientation of B in terms of the polar 
and azimuthal angles, 8 and 0, measured with respect to the 
axes of the magnetic complex. If a value for 0 is not given, 
then the Hamiltonian is axially symmetric. For materials with 
non-equivalent sites, 8 and 0 are restricted to values for which 
the sites are equivalent. The fifth row of the table gives the 
magnitude of B. The fields do not exceed 14 kG for the oper- 
ating points shown. 

Rows 6,7, and 8 list the signal and pump frequencies with 
the corresponding transition levels shown in parentheses. The 
signal frequency is always 32.0 GHz. The pump frequencies 
vary roughly between 50 and 90 GHz. Larger pump frequencies 
will in general yield larger values of I opt and l equaV On the 
other hand, copper and dielectric losses increase at higher fre- 
quencies, with the result that heating of the maser structure 
may restrict the use of high pump power levels at high fre- 
quencies for some materials. Note also that pump frequencies 
within the same waveguide band simplify engineering issues. 

Row 9 gives the value of AfJAB evaluated near the operat- 
ing point. For a maser tunable over a wide range, AfJAB 
should be of the same sign and of similar magnitude for the 
maser material and the isolator material. Rows 10 and 11 give 
values of Af p jAf s (actually, Af p /AB • AB/Af s ) for both 
pumps. This parameter is indicative of the pump bandwidth 
required for a given signal bandwidth, so it is preferable for 
| Af p /Af s | to be as small as possible. For most of the operating 
points Af p /Af s » 2, but several have values <1. Note that 
Af/AB is evaluated as a simple two-point difference, with the 
second point arbitrarily located 200 G from the operating 
point. For strongly curved energy levels, these values may not 
be accurate across the desired band. 

Rows 12, 13, and 14 give a 2 for the signal and pump transi- 
tions for optimum elliptieally polarized fields according to 
Eq. (10). Recall from Eqs. (14) and (16) that G dB and the 
figure-of-merit are proportional to a 2 , so as large a value of a 2 
as possible is desired. In general, a 2 is a factor of 2- to 3-times 
larger for the 6-level spin systems. Similarly, a large value of 
a 2 is preferred, since the pump power required for saturation 


is inversely proportional to a 2 . According to [29] , the pump 
power required for saturation will satisfy 


P cc 

pump 


L • A/„ 


a 2 t 
P P 


(17) 


where r p is the effective pump relaxation time;r p is not iden- 
tical to the measured pump relaxation time. The values of a 2 
in Table 1 span nearly two orders of magnitude. 


Rows 15 and 16 show the inversion ratios for equal and 
optimum relaxation times. The values of I equai and I opt are 
similar for the various operating points, except when only one 
pump is employed. 

Finally, row 17 gives the material^ figure-of-merit, com- 
puted in units of MHz, as 


A^ = 5.6/ op ,-q, 2 

Q no. of levels 

m 


(18) 


This follows from Eq. (13) evaluated at/= 32 GHz, T= 4.2 K, 
7] = 0.5, and A" = 2.35 X 10 19 spins/cm 3 . Since the true inver- 
sion ratio for a given operating point may be as much as a 
factor of 3 or more smaller than I opv a detailed comparison 
of figure-of-merits could be misleading. 


Table 1 is by no means complete in the sense that one may 
confidently select the best maser material from it. Two very 
important parameters are missing: the actual inversion ratio 
and the pump power, P pump , required to maintain that inver- 
sion ratio. At present, both of these parameters must be 
measured. 


Since both / and P pump depend critically on relaxation 
times, a table of relaxation times, located in the Appendix and 
labeled Table A, was compiled from data found in the litera- 
ture. Ionic concentration, frequency, orientation, and transi- 
tion information is included. Because these parameters do not 
coincide with pump operating-points of interest to us, and 
because of the dependence of relaxation times on measure- 
ment technique and crystal growth procedures [30] , the data 
in Table A could easily be an order of magnitude or more 
different from what would be measured for the materials in 
Table 1 . Hence, the relaxation times in Table A are not used in 
any calculations in this work, even though they are the best 
values available to us at the present time. 

Before discussing the many materials in Table 1 , consider 
the operating point in which ruby is presently being used at 
32 GHz (first column of table, Ruby No. 1). The ruby is 
oriented at the double-pump angle (8 = 54.7 deg) and pumped 
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in the push-pull mode, so the pump frequencies are equal. 
Scanning down the column, two potential problems can be 
seen with this operating point. First, the pump bandwidths 
are nearly twice the signal bandwidths, so if 500 MHz of 
signal bandwidth is desired at 32 GHz, 1-GHz bandwidth must 
be pumped at 66 GHz. To pump such a large bandwidth, the 
pumps must be swept across the band, effectively reducing the 
pump power at a given frequency. How detrimental this is 
depends on the relaxation times of the pump transitions. 

The second problem with ruby at this operating point is 
the weak pump transitions: o 2 x - 0.05 and o 2 2 = 0.04 com- 
pared to o 2 = 1.92. For this reason, high levels of pump power 
are used in the 32-GHz RWM, although the pumped levels are 
still not saturated. Note that a small o 2 does not preclude 
good maser performance, as demonstrated by the 18- to 
26-GHz maser of Moore and Clauss [2], [31] for which 

°h = and ° 2 P 2 = 0-11- 

The most significant problem with ruby, the low inversion 
ratio, is not indicated by the table. Measured values of / for 
the case of saturated pump transitions have been approxi- 
mately 1.1 ([3], and J. Shell, private communication). (The 
similarity to I equaX = 1.1 does not necessarily mean that the 
relaxation times are equal.) In the 18- to 26-GHz range, mea- 
sured values of / have been in the range of 1.6 to 1.8 [2] , 

[3] > [31] * 

Finally, from Table A it can be seen that ruby has long 
relaxation times compared to the other materials. Hence, even 
though o 2 is small, the denominator of Eq. (17) remains large 
enough for ruby to require large but manageable pump power. 

Consider several other operating points in Table 1. Since 
ruby has worked so well in the past, ruby at another orienta- 
tion is an obvious candidate for a maser material. The second 
column of the table, Ruby No. 2, shows ruby at 9 = 90 deg 
and with a push-push pumping scheme. Even though a 2 and 
o p j are weaker and the values for I equal and I opt are less than 
for Ruby No. 1, if the actual inversion ratio is >1.5, Ruby 
No. 2 could yield a higher gain -bandwidth product. Some 
investigators [4] claim Ruby No. 2 to be supejior to Ruby 
No. 1 at millimeter wavelengths because of a higher inversion 
ratio and less critical orientation (less spreading of pump 
power due to c-axis wander). 

The sapphire host has many desirable properties, so Fe- 
doped sapphire is a logical choice. For Sapphire No. 1, o 2 is 

2- to 5-times stronger than that of Ruby No. 1 and o ] is 2- to 

3- times that of Ruby No. 1. However, according to Table A, 
the pump relaxation times may be an order of magnitude 
shorter, implying that Sapphire No. 1 could require several 
times the pump power of Ruby No. 1. Other investigators 


[28] suggest that the relaxation times of Fe-doped sapphire 
are similar to those of ruby. If this is true, then Sapphire No. 1 
could require several times less pump power than Ruby No. 1 . 
Measurements of the relaxation times and pump power re- 
quired must be made to determine which scenario is correct. 

Emerald has some similarity to ruby, having the same spin 
Hamiltonian and potentially long relaxation times. If the 
inversion ratio for Emerald No. 1 is >2, then this operating 
point would be very attractive. A problem with emerald is the 
difficulty of its growth, which may not allow the high degree 
of crystal perfection necessary. 

Zinc tungstate has a complicated H s , large ZFS values, and 
may have short relaxation times, making it quite different 
from ruby. Cr-doped zinc tungstate has several promising oper- 
ating points. In particular, ZnW0 4 No. 1 is attractive, assum- 
ing I ^ l opr Fe-doped zinc tungstate exhibits a large number 
of excellent operating points. For ZnW0 4 No. 3, a 2 is 2- to 
3-times greater than that of Ruby No. 1 and o 2 is 100 times 
that of Ruby No. 1. This large value of o 2 raises the question 
of whether it is preferable to have large a 2 and small r p 
(ZnW0 4 No. 3) or small o 2 and large r p (Ruby No. 1). Assum- 
ing the product o 2 * r p is constant, the pump power require- 
ments will be similar, but in the former case more energy 
would be transferred to the lattice. This could raise the tem- 
perature of the maser material, thereby decreasing the gain; 
however, this possibility has not been considered in detail. 
Harmonic cross relaxation may be a problem for several of the 
better operating points for ZnW0 4 :Fe. One similarity zinc 
tungstate has with ruby is that it can be grown by the Czoch- 
ralski method. 

The rutiles appear promising, but the large, anisotropic, 
temperature-dependent dielectric constant of rutile makes it 
unattractive from an engineering standpoint. 

One can easily see from Table 1 that many of the other 
materials analyzed may make excellent maser materials, but 
the lack of information on inversion ratios, pump power 
requirements, relaxation times, etc., makes them difficult to 
evaluate. 

Another possibility, not addressed in Table 1, is to use 
standard ruby doped with a fast-relaxing impurity. This addi- 
tional impurity may be added in the melt or created in the 
finished ruby by exposure to X-rays (so-called orange-ruby 
[9]). A properly chosen impurity can shorten certain relax- 
ation times, which, by making the times more optimal, can 
increase the inversion ratio. However, the impurity would not 
alter the ZFS, so the pump transitions would still be weak. 

We hope to eventually make measurements of inversion 
ratios, pump power requirements, and relaxation times at 
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32 GHz and around 60 GHz on several of the materials in 
Table 1. 

Other materials we would like to analyze but for which 
we do not have the spin Hamiltonians are spinekFe and 
chrysoberyhCr, Fe. 

V. Conclusions and Future Work 

Any of the materials analyzed in this work may yield better 
maser performance than does ruby at 32 GHz at the double- 
pump angle. However, several key parameters related to pump 
power requirements may eliminate some or all of these mate- 
rials. Based on results from the analysis of the spin Hamil- 
tonians and on scanty (and unreliable) relaxation time data, 
several materials show particular promise (e.g., Fe-doped zinc 
tungstate). 


To complete the materials evaluation, it will be necessary to 
measure the inversion ratio and pump power required for 
saturation for each operating-point of interest. Barring cross- 
relaxation and other concentration-dependent effects, know- 
ledge of the relaxation-times would be sufficient to calculate 
both / and P pump . However, the subtlety of measuring 
relaxation-times will most likely require that / and P pump be 
measured. 


Better understanding of the low inversion ratio of ruby is 
needed. By accounting for the spin-phonon interaction, one 
can calculate the relaxation rates of the transitions for low 
spin concentrations [29] . With these relaxation rates, the 
inversion ratios and pump power requirements can be calcu- 
lated for each operating-point of interest and for various 
physical temperatures. 
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Table 1. Promising 32-GHz operating points for ruby and other materials. The materials are 

arranged in order of increasing ZFS. 





Material 



Operating 













Point 

Ruby 

Ruby 

Ruby 

Sapphire 

Sapphire 

YAG 


No. 1 

No. 2 

No. 3 

No. 1 

No. 2 

No. 1 

Ion 

Cr 

Cr 

Cr 

Fe 

Fe 

Cr 

ZFS, GHz 

11.4 

11.4 

11.4 

12.1, 19,1 

12.1, 19.1 

15.7 

No. of ionic sites 

1 

1 

1 

2 

2 

1 

0,0, deg 

54.74 

90 

90 

90,45 

60,30 

54.74 

B, kG 

11.81 

13.50 

11.20 

9.50 

12.44 

13.27 

4.GHz 

32.0 (32) 

32.0(21) 

32.0 (32) 

32.0 (32) 

32.0 (54) 

32.0 (32) 

fp\ 

66.2 (13) 

70.3(13) 

57.6 (13) 

68.6 (13) 

77.3(13) 

76.2(13) 

fp2 

66.2 (24) 

43.3 (34) 

68.9 (24) 

59.9 (24) 

66.1 (35) 

76.2 (24) 

AfjAB, MHz/G 

2.9 

2.8 

2.7 

2.8 

3.0 

2.9 

Afpi/Af, 

1.9 

2.0 

2.0 

2.0 

1.8 

1.8 


1.9 

1.0 

2.0 

2.0 

2.0 

1.8 


1.92 

1.51 

1.97 

6.87 

5.64 

1.62 

°l l 

0.05 

0.02 

0.03 

0.08 

0.09 

0.13 

°P2 

0.04 

1.51 

0.02 

0.17 

0.28 

0.31 

I opt 

3.1 

2.2 

2.8 

3.0 

4.8 

3.7 

1 equal 

1.1 

0.7 

0.9 

1.0 

2.0 

1.4 

A4/fi m ,MHz 

8.3 

4.7 

7.7 

19.2 

25.3 

8.4 




Material 



Operating 













Point 

YAG 

YGG 

YGG 

Spinel 

Andalusite 

Andalusite 


No. 2 

No. 1 

No. 2 

No. 1 

No. 1 

No. 2 

Ion 

Cr 

Cr 

Cr 

Cr 

Cr 

Cr 

ZFS, GHz 

15.7 

20.9 

20.9 

29.7 

32.0 

32.0 

No. of ionic sites 

1 

1 

1 

4 

2 

2 

0,0, deg 

70 

54.74 

70 

54.74 

55,0 

70, 0 

B y kG 

13.40 

13.96 

13.35 

13.15 

12.94 

13.00 

fs' GHz 

32.0 (43) 

32.0 (32) 

32.0 (43) 

32.0 (32) 

32.0 (32) 

32.0 (43) 

fpl 

49.6 (12) 

81.8(13) 

53.3 (12) 

75.3 (13) 

74.2(13) 

47.0(12) 

fp 2 

69.7 (24) 

81.8 (24) 

70.7 (24) 

75.3 (24) 

75.8 (24) 

69.4 (24) 

Af s /AB, MHz/G 

2.4 

2.7 

2.2 

2.9 

2.8 

2.4 

A/p i /A/, 

1.2 

1.9 

1.3 

1.8 

1.9 

1.1 

Af p2 /Af s 

2.1 

1.9 

2.1 

1.8 

1.9 

2.1 

•? 

1.25 

1.50 

1.01 

1.66 

1.63 

1.27 

°|l 

1.52 

0.16 

1.54 

0.28 

0.15 

1.53 

°P2 

0.31 

0.47 

0.55 

0.13 

0.31 

0.31 

1 opt 

3.6 

4.0 

3.9 

3.6 

3.6 

3.5 

I equal 

1.5 

1.6 

1.6 

1.4 

1.4 

1.4 

Af L !Q m yUHz 

6.3 

8.4 

5.5 

8.4 

8.2 

6.2 
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Table 1 (contd) 


Operating 

Point 



Material 



Rutile 
No. 1 

Rutile 
No. 2 

Rutile 
No. 3 

Rutile 
No. 4 

ZnW0 4 
No. 1 

ZnW0 4 
No. 2 

Ion 

Cr 

Cr 

Fe 

Fe 

Cr 

Cr 

ZFS, GHz 

43.3 

43.3 

43.3,81.3 

43.3,81.3 

51.6 

51.6 

No. of ionic sites 

2 

2 

2 

2 

1 

1 

<9, 0 , deg 

45,0 

54.74,45 

52.55,40 

71.12,70 

40,90 

50,90 

B t kG 

12.78 

14.06 

9.35 

11.71 

9.78 

13.10 

GHz 

32.0 (43) 

32.0 (32) 

32.0 (32) 

32.0 (43) 

32.0(21) 

32.0 (32) 

fpl 

56.0(12) 

82.6 (13) 

78.9 (13) 

71.8 (12) 

54.1 (13) 

75.1 (13) 

fp2 

65.4 (24) 

82.6 (24) 

73.3 (24) 

81.4 (24) 

54.9 (34) 

88.4 (24) 

Af s /AB, MHz/G 

2.0 

2.6 

2.8 

3.7 

2.3 

1.9 

AApi/AT, 

1.5 

1.9 

2.1 

0.9 

1.1 

2.2 

&fp2^fs 

2.2 

1.9 

2.1 

2.1 

1.8 

2.6 

° 2 S 

1.57 

1.37 

3.80 

3.57 

1.68 

1.37 


1.18 

0.48 

0.67 

3.31 

0.53 

0.31 

°P2 

0.37 

0.30 

1.55 

1.01 

0.57 

0.61 

I opt 

3.7 

4.1 

3.7 

5.7 

2.1 

3.8 

I equal 

1.4 

1.6 

1.4 

2.6 

0.5 

1.5 

Af L,IQm’ MHz 

8.1 

7.9 

13.1 

19.0 

4.9 

7.3 




Material 



Operating 













Point 

Emerald 

Emerald 

ZnW0 4 

ZnW0 4 

Y 2°3 

Andalusite 


No. 1 

No. 2 

No. 3 

No. 4 

No. 1 

No. 3 

Ion 

Cr 

Cr 

Fe 

Fe 

Cr 

Fe 

ZFS, GHz 

53.5 

53.5 

61.0,76.9 

61.0,76.9 

72.7 

112.6,225.2 

No. of ionic sites 

1 

1 

1 

1 

4 

2 

6 , 0 , deg 

40,0 

54.74 

90,45 

90,45 

40 

30,0 

B, kG 

7.95 

14.04 

8.59 

5.28 

9.14 

6.63 

/..GHz 

32.0 (43) 

32.0 (32) 

32.0 (54) 

32.0 (54) 

32.0 (21) 

32.0 (21) 

fp 1 

48.4 (12) 

86.1 (13) 

66.0(13) 

67.9(13) 

69.1 (13) 

107.3 (13) 

fp2 

50.8 (24) 

86.1 (24) 

67.4 (35) 

57.6 (35) 

88.3 (24) 

- 

Af s /AB, MHz/G 

3.5 

2.2 

2.1 

-2.7 

-3.32 

4.78 


1.5 

2.2 

-0.1 

0.4 

-0.16 

-0.07 

Af p2 lAf s 

0.3 

2.2 

2.1 

-0.5 

-0.72 

- 

°S 

1.54 

1.38 

5.52 

3.24 

0.76 

3.55 

°P 1 

0.13 

0.66 

3.01 

3.17 

0.73 

2.32 

°h 

0.59 

0.19 

2.91 

3.51 

1.31 

- 

I opt 

2.6 

4.3 

4.3 

3.8 

3.6 

2.0 

I equal 

0.8 

1.7 

1.8 

1.5 

1.4 

0.5 

A4/e m ,MHz 

5.6 

8.3 

22.2 

11.5 

3.8 

6.6 
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Appendix 

Table of Relaxation Times 


This Appendix contains a table of relaxation times (7^) for 
various materials at 4.2 K. The paramagnetic ion concentra- 
tion, transition frequency, transition, and orientation are also 
shown. In addition, the literature references are included. Note 
that the table does not mention either the crystal growth 


process or the relaxation time measurement technique, both of 
which may significantly alter the reported relaxations times. 
Hence, these relaxation times could easily be an order of mag- 
nitude or more different from what one might measure, and 
should not be used in calculations unless verified. 


Table A-1 . Measured relaxation times, T r for various materials at physical temperature T = 4.2° K 


Material 

Ionic 

concentration, 
atomic percent 

Frequency, GHz 

Orientation, deg 
(9,0) 

Transition 

Tj, msec 

Reference 

Ruby 

0.03 

34.6 

90 

3-4 

21 

[32] 





2-3 

16 






1-2 

22 






2-4 

54 






1-3 

56 



0.013 

35 

90 

2-3 

15.5 

[33] 


0.052 



2-3 

17.5 


Sapphire: Fe 

0.03 

34.6 

90,0 

4-5 

1.8 

[34] 





3-4 

1.6 






2-3 

1.5 






1-2 

2.0 


Emerald 

4.9 X 10 19 ions/cm 3 

9.3 

0 

3-4 

9 

[35] 




90 

1-2 

8 






3-4 

11 


ZnW0 4 : Cr 

0.005-0.3 

9.2 

90, 90 

1-2 

*1.5 

[36] 

0.018-0.72 

33 

90,90 

1-2 

*0.5 



0.08 

X-band 

? 

5-6 

~0.3 

a 





1-2 

-1 


Rutile: Cr 

0.07 

34.6 

90, 0 b 

1-4 } SitCA 

4.5 

2.5 

[34] 





il } siteB 

2.5 

2.1 


Rutile: Fe 

0.01-0.02 

9.4 

0,0 

1-2 

*2 

[37] 





3-4 

*2 


a J. Orton, private communication with K. Standley and R. Vaughan. 
b For this orientation the ionic sites are inequivalent. 
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